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ABSTRACT 

We present spatially-resolved X-ray spectroscopy of the luminous lensing cluster Abell 
2390, using observations made with the Chandra observatory. The temperature of the 
X-ray gas rises with increasing radius within the central ~ 200 kpc of the cluster, and 
then remains approximately isothermal, with kT = 11.5^'g keV, out to the limits of 
the observations at r ~ 1.0 Mpc. The total mass profile determined from the Chandra 
data has a form in good agreement with the predictions from numerical simulations. 
Using the parameterization of Navarro, Frenk & White (1997), we measure a scale 
radius r s ~ 0.9 Mpc and a concentration parameter c ~ 4. The best-fit X-ray mass 
model is in good agreement with independent gravitational lensing results and optical 
measurements of the galaxy velocity dispersion in the cluster. The X-ray gas-to-total- 
mass ratio rises with increasing radius with / gas = 21 ± 10 per cent at r = 0.9 Mpc. 
The azimuthally-averaged 0.3 — 7.0 keV surface brightness profile exhibits a small core 
radius and a clear 'break' at r ~ 500 kpc, where the slope changes from Sx 9S r~ 1,5 
to Sx SS r~ 3 6 . The data for the central region of the cluster indicate the presence of 
a cooling flow with a mass deposition rate of 200 — 300 M Q yr _1 and an age of 2 — 3 
Gyr. 
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1 INTRODUCTION 

Accurate measurements of the masses of clusters of galaxies 
are of profound importance to cosmological studies. Origi- 
nally, most measurements of cluster masses were based on 
optical studies of their galaxy dynamics, wherein the mo- 
tions of individual galaxies were used to trace the cluster 
potentials. Although such studies were shown to be sensi- 
tive to systematic uncertainties due to velocity anisotropies, 
substructure and projection effects (e.g. Lucey 1983; Frenk 
et al. 1990; van Haarlem, Frenk & White 1997), more recent 
work based on large galaxy samples and employing careful 
selection techniques, has lead to significant progress (e.g. 
Carlberg et al. 1996; den Hartog & Katgert 1996; Fadda et 
al. 1996; Mazure et al. 1996; Borgani et al. 1999; Geller, 
Diaferio & Kurtz 1999; Koranyi & Geller 2000). 

Recent years have also seen the development of two fur- 
ther techniques for measuring the masses of clusters, based 
on X-ray observations and studies of gravitational lensing 
by clusters, respectively. X-ray mass measurements use the 
assumption that the X-ray emitting gas which pervades clus- 
ters is in hydrostatic equilibrium; the total mass distribution 
is determined once the radial distributions of the X-ray gas 
density and temperature are known (see e.g. Sarazin 1988). 
Since the X-ray emissivity is proportional to the square of 



the gas density, and the relaxation timescale for the X-ray 
gas is relatively short (of the order of a few sound crossing 
times), the X-ray method is relatively free from the projec- 
tion and substructure effects which hamper the aforemen- 
tioned optical studies. 

In contrast to the X-ray and optical dynamical tech- 
niques, gravitational lensing offers a method for measuring 
the projected masses through clusters that is essentially free 
from assumptions about the dynamical state of the grav- 
itating matter (see e.g. Fort & Mellier 1994; Bartelmann 
& Schneider 1999; Mellier 1999 for reviews). The primary 
observational challenges of requiring deep exposures, ex- 
cellent seeing conditions, wide field imaging and accurate 
point spread function models have now been mostly over- 
come with improved instrumentation (e.g. Bacon et al. 2000 
and references therein), although the recovery of the three- 
dimensional mass distributions in clusters can be compli- 
cated by projection effects and uncertainties in the redshift 
distributions of the lensed sources. 

Clearly, the best approach when attempting to reliably 
measure the masses of clusters is to combine these three 
methods. The first combined X-ray and lensing studies of 
galaxy clusters (Miralda-Escude & Babul 1995) suggested 
that strong lensing masses, measured within r < 200 kpc 
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of the cluster centre, might typically overestimate X-ray- 
determined masses by a factor of ~ 2 — 3. This sparked 
debate into the possible effects of oblate/prolate cluster ge- 
ometries, projection effects, complex temperature structures 
and pressure support from bulk and/or turbulent motions 
and magnetic fields in the X-ray gas (e.g. Miralda-Escude 
& Babul 1995; Loeb & Mao 1994; Waxman & Miralda- 
Escude 1995; Kneib et al. 1995; Bartelmann & Steinmetz 
1996). Later work (Allen, Fabian & Kneib 1996; Allen 1998; 
Bohringer et al. 1998; Wu 2000) highlighted a clear differ- 
ence between the results obtained for cooling-flow (hereafter 
CF) and non-cooling flow (NCF) clusters. For CF clusters, 
the X-ray and strong lensing mass measurements generally 
show good agreement, once the effects of the cooling flows 
are accounted for in the X-ray analysis. For NCF systems, 
however, the masses inferred from the strong lensing data 
invariably exceed the X-ray values, determined under the 
hydrostatic assumption, by a factor of 2 — 4. 

The origin of the different results obtained for CF and 
NCF clusters is thought to lie in the different dynamical 
states of these systems: whereas X-ray and optical imaging 
and detailed lensing analyses of CF clusters show them to 
be relatively regular and dynamically-relaxed systems, NCF 
clusters generally appear to be undergoing major subcluster 
merger events (e.g. Edge, Stewart & Fabian 1992; Buote & 
Tsai 1996; Kneib et al. 1995; Smail et al. 1995, 1997; Squires 
et al. 1997). Significant offsets between the X-ray and lensing 
centroids are observed in NCF clusters, demonstrating a loss 
of hydrostatic equilibrium in their central regions. The X- 
ray core radii for NCF systems also appear to have been 
inflated by the dynamical activity, in agreement with the 
predictions from numerical simulations (e.g. Roetigger et al. 
1996); this inflation of the X-ray core radii can account for 
the bulk of the X-ray/strong lensing mass discrepancy in 
most NCF systems (Allen 1998). On larger spatial scales 
(r > 0.5 Mpc), comparisons between weak lensing, X-ray 
and optical dynamical mass measurements generally provide 
a consistent picture, with excellent agreement between the 
results obtained for both CF and NCF clusters (e.g. Squires 
et al. 1996; Smail et al. 1997; Wu & Fang 1997; Allen 1998; 
Lewis et al. 1999). This suggests that the loss of hydrostatic 
equilibrium in NCF clusters is primarily restricted to their 
inner regions. 

The most significant uncertainty associated with the 
X-ray mass measurements in previous joint X-ray/lensing 
studies has been the absence of any direct measurements 
of the X-ray temperature profiles in the clusters, which im- 
pacts directly on the mass measurements through the hy- 
drostatic equation. Although limited spatially-resolved spec- 
troscopy for bright, nearby clusters was possible using ASCA 
and Beppo-SAX observations (e.g. Markevitch et al. 1998; 
Kikuchi et al. 1999; de Grandi & Molendi 1999; White 2000; 
Irwin & Bregman 2000), for the more distant lensing clus- 
ters, typically observed at redshifts z > 0.2, only a single, 
integrated cluster spectrum was normally available. In most 
cases, only a mean emission-weighted X-ray temperature 
was therefore determined, although more sophisticated stud- 
ies also accounted for the effects of cooling flows on the X-ray 
data (e.g. Allen 1998; Bohringer et al. 1998). For the mass 
analyses, it was then necessarily assumed that the mass- 
weighted temperature profile followed some particular form, 
and usually that it remained approximately isothermal with 



radius. However, the validity of this assumption remains un- 
certain, especially within the strong lensing regime. 

The launch of the Chandra Observatory (Weisskopf et 
al. 2000) in 1999 July provides the first opportunity for de- 
tailed, spatially-resolved X-ray spectroscopy of clusters of 
galaxies at moderate redshifts. The Advanced CCD Imag- 
ing Spectrometer (ACIS) on Chandra permits the first di- 
rect, simultaneous measurements of the X-ray temperature 
and density profiles and, via the hydrostatic assumption, 
the mass distribution in luminous lensing clusters, spanning 
both the weak and strong lensing regimes. In this paper we 
present the first results from Chandra observations of the 
massive CF lensing cluster Abell 2390, which has been the 
subject of several previous combined X-ray /optical/lensing 
studies (e.g. Piere et al. 1996, Squires et al. 1996, Allen 1998, 
Bohringer et al. 1998, Lewis et al. 1999). We present detailed 
results on the mass distribution in the cluster determined 
from the Chandra data and compare our results with those 
from detailed strong and weak lensing analyses and optical 
dynamical studies. We also re-examine the properties of the 
cooling flow in the cluster and discuss their relation to the 
dynamical history of the system. 

The cosmological parameters _Ho=50 kms -1 Mpc - , 
Q = 1 and A = are assumed throughout. At the red- 
shift of Abell 2390 (z — 0.2301) an angular scale of 1 arcsec 
corresponds to a physical length of 4.652 kpc. 



2 OBSERVATIONS 

The Chandra observations of Abell 2390 were carried out 
using the ACIS on 1999 November 7. The target was ob- 
served in the back-illuminated CCD detectors, close to the 
nominal aim point for the ACIS S3 detector. (The source 
was positioned near the centre of node-1 on the chip). The 
focal plane temperature at the time of the observations was 
-HOC. 

We have used the CI AO software (version 1.1.3) and 
the level-2 events file provided by the standard Chandra 
pipeline processing for our analysis. The light curve for the 
observation was of high quality with no strong background 
flaring. Only those X-ray events with grade classifications of 
0,2,3,4 and 6 were included in our final cleaned data set, for 
which the net exposure time was 9.13ks. 



3 X-RAY IMAGING ANALYSIS 
3.1 X-ray morphology 

The raw 0.3 — 7.0 keV image of the central 6x6 arcmin 2 
(1.8 x 1.8 Mpc 2 ) region of Abell 2390 is shown in Fig. 
h|(a). The pixel size is 1.97 x 1.97 arcsec 2 , corresponding 
to 4 x 4 raw detector pixels. Fig. [l](b) shows an adaptively 
smoothed version of the same image, using the smoothing 
algorithm of Ebeling, White & Rangarajan (2000). The po- 
sition of the peak of the X-ray emission from the cluster, 
21h53m36.77 +17d41m42.8s (J2000.), is in excellent agree- 
ment with the optical centroid for the dominant cluster 
galaxy of 21h53m36.76 +17d41m42.9s (J2000.; Pierre et al. 
1996). The X-ray image is elongated along an approximately 
northwest-southeast direction, in a similar manner to the op- 
tical isophotes of the dominant galaxy (Pierre et al. 1996). 



© 0000 RAS, MNRAS 000, 000-000 



Chandra measurements of the distribution of mass in Abell 2390 3 



Some substructure is apparent in the Chandra image. 
In the central region, a relatively bright ridge of enhanced 
emission extends ~ 3 — 4 arcsec to the northwest of the X- 
ray peak (Fig. ^) . Enhanced emission in the same direction is 
also observed in the optical blue continuum and optical/UV 
line emission from the dominant cluster galaxy (Lemonon 
et al. 1998; Edge et al. 1999; Hutchings & Balogh 2000). A 
correspondence between the brightest X-ray emission with 
excess blue continuum and optical/UV line emission (prob- 
ably associated with the formation of young, massive stars) 
is also observed in the nearby CF clusters Hydra-A (McNa- 
mara et al. 2000) and Abell 1795 (Fabian et al. 2000b). A 
second, fainter region of enhanced X-ray emission extends 
~ 6 arcsec to the south-southeast. 

On medium (r ~ 40 — 80 arcsec) scales, the X-ray emis- 
sion is extended to the west-northwest, in a similar manner 
to the optical luminosity distribution and lensing mass mod- 
els of Piere et al. (1996; these authors also note the presence 
of an X-ray extension in the same direction using ROSAT 
High Resolution Imager data). On large scales (> 2 arcmin) 
the X-ray emission is extended towards the east. The pres- 
ence of such substructure in the Chandra image suggests 
that the cluster has not fully relaxed following it's most re- 
cent merger activity. However, the agreement between the 
X-ray, optical and gravitational lensing mass measurements 
discussed in Section 5.2 argues that overall the assumption 
of hydrostatic equilibrium in the cluster is a reasonable one. 

Three point sources are detected to the west of the X- 
ray peak, in the central regions of the cluster, at positions 
21h53m33.19 +17d42m08.4s, 21h53m33.74 +17d41ml3.3s, 
and 21h53m34.02 +17d42m39.9s. The X-ray and sub-mm 
properties of these sources are discussed by Fabian et al. 
(2000a). Four more point sources are also visible at larger 
radii in Fig. |l|(b). 



3.2 The surface brightness profile 

The azimuthally-averaged, 0.3 — 7.0 keV X-ray surface 
brightness profile for Abell 2390 is shown in Fig. ^. The 
profile has been flat-fielded and background subtracted us- 
ing a rectangular background region of size 400 x 50 arcsec 2 , 
located ~ 5 arcmin from the cluster centre. All obvious point 
sources were excluded from the analysis. 

The X-ray emission from the cluster extends beyond the 
5 arcmin (~ 1.4 Mpc) radius associated with our on-chip 
background region (see also Pierre et al. 1996; Bohringer et 
al. 1998). However, beyond this radius background counts 
dominate the flux in the ACIS-S3 detector. Fig. ^| shows the 
data for the central 900 kpc (193 arcsec), for which system- 
atic errors associated with the background subtraction and 
flat fielding are negligible. The bin-size in Fig. ^(a) is 2 de- 
tector pixels (0.984 arcsec) . Fig. ^(b) shows the data for the 
outer regions of the cluster with a larger binsize of 8 detector 
pixels (3.94 arcsec). 

Within a radius of 500 kpc (107 arcsec), the X-ray sur- 
face brightness profile can be parameterized (x 2 = 119 for 
106 degrees of freedom) by a standard /3-model (e.g. Jones & 
Forman 1984) of the form S(r) = 5(0) [(1 + r/r c ) 2 ] 1/2_3 ' 3 ) 
with a core radius r c = 28. 7± 1.3 kpc and a slope parameter 
/3 = 0.412 ± 0.003 (Iff errors; Ax 2 = 1.0). On larger scales, 
however, the /3-model does not provide an acceptable fit: ex- 



amining the data for the central 900 kpc radius, we obtain 
X 2 = 909 for 194 degrees of freedom, with best-fit parameter 
values of r c = 55 kpc and /3 = 0.48). Ignoring the central 
r ~ 80 kpc (18 arcsec) region, associated with the possi- 
ble cooling flow (Section 6), the fit is improved (x 2 = 485 
for 177 degrees of freedom, with best-fit parameter values 
r c = 160kpc and ft = 0.58) although is still formally unac- 
ceptable. 

The main reason for the poor fit obtained with the j3- 
model at larger radii is the presence of a 'break' in the sur- 
face brightness profile at r ~ 500 kpc (Fig. ^b). This break 
is not obviously due to substructure in any particular di- 
rection in the cluster: Fig. ^ shows the surface brightness 
profile measured in the four quadrants covering position an- 
gles 45 - 135, 135 - 225, 225 - 315 and 315 - 45 degrees. 
We see that the profile appears remarkably similar in three 
of the four directions, although the emission is slightly more 
extended towards the east (as is also evident in the images 
presented in Fig. [j]). 

A good fit to the surface brightness profile beyond the 
central cooling region can be obtained using a simple broken 
power- law model. Fitting the data from 80 — 900 kpc, we ob- 
tain x 2 = 173 for 176 degrees of freedom, with a break at a 
radius of 49lly 6 kpc, and slopes in the regions internal and 
external to the break radius of — 1.46±0.02 and — 3.60±0.17, 
respectively. Interestingly, these slopes are similar to the val- 
ues expected at small and large radii for the dark matter in 
a Navarro, Frenk & White (1997; hereafter NFW) potential 
in which poc l/((r/r s )(l+r/r s ) 2 ). However, isothermal gas 
in an NFW-like potential should not exhibit a sharp break 
at the scale radius r s , but rather a slow rollover (although 
individual clusters in the simulations presented by Thomas 
et al. 2000 do exhibit sharp breaks in their dark matter dis- 
tributions). Fitting the surface brightness profile external 
to the cooling flow (80 — 900 kpc) with the prescription for 
isothermal gas in an NFW potential described by Ettori & 
Fabian (1999), we measure r s = 610±f kpc (formal Iff er- 
rors), with x 2 = 420 for 177 degrees of freedom. Thus, the 
fit with the NFW mass model assuming strict isothermality 
in the X-ray gas is formally unacceptable, although provides 
a better description of the data in the 80 — 900 kpc region 
than the /3-model. As discussed in Section 4, the Chandra 
data show that the X-ray gas in Abell 2390 is not isothermal 
and that the temperature rises with increasing radius within 
the central r ~ 200 kpc. In Section 5 we show that an NFW 
mass model can provide a good description of the Chandra 
data, once the assumption of isothermality is relaxed. 



4 SPATIALLY-RESOLVED SPECTROSCOPY 
4.1 Method of analysis 

For our spectral analysis, we divided the cluster into annu- 
lar regions, as detailed in Table 1. A spectrum was extracted 
from each region in 1024 Pulse Height Analyser (PHA) chan- 
nels. The spectra were re-grouped to contain a minimum of 
20 counts per PHA channel, thereby allowing x 2 statistics 
to be used. (For the two outer annuli, a larger grouping of 
40 counts per PHA channel was used, due to the increased 
background contribution). Background spectra, appropriate 
for the regions studied, were extracted from the ACIS-S3 
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blank-field data sets available from the Chandra X-ray Cen- 
ter. All obvious point sources were masked out and excluded 
from the analysis. Separate photon-weighted response ma- 
trices and effective area files were constructed for each region 
using the calibration and response files appropriate for the 
focal plane temperature, available from the Chandra X-ray 
Centre. 

Two separate energy ranges were examined. Firstly, a 
conservative 0.5 — 7.0 keV band was used over which the 
calibration of the back-illuminated CCD detectors is cur- 
rently best understood. Secondly, for the central 100 kpc 
region, where a cooling flow is thought to exist (Section 6), 
we have also examined a more extended 0.3 — 7.0 keV energy 
range, which provides extra constraints on the presence of 
cool emission components and/or intrinsic absorption in the 
cluster. 

4.2 The spectral models 

The analysis of the spectral data has been carried out us- 
ing the XSPEC software package (version 11.01; Arnaud 
1996). The spectra were modeled using the plasma emis- 
sion code of Kaastra & Mewe (1993; incorporating the Fe 
L calculations of Liedhal, Osterheld & Goldstein 1995) and 
the photoelectric absorption models of Balucinska-Church & 
McCammon (1992). We first examined each annular spec- 
trum using a simple, single-temperature model with the ab- 
sorbing column density fixed at the nominal Galactic value 
(Nn = 6.8 x 10 20 atom cm" 2 ; Dickey & Lockman 1990). This 
model is hereafter referred to as model A. The free parame- 
ters in model A were the temperature (kT) and metallicity 
(Z) of the plasma (measured relative to the solar photo- 
spheric values of Anders & Grevesse 1989, with the various 
elements assumed to be present in their solar ratios) and the 
emission measure (K). We also examined a second single- 
temperature model (model B) which was identical to model 
A but with the absorbing column density (Nn) also included 
as a free parameter in the fits. 

The image deprojection and X-ray colour profile analy- 
ses discussed in Section 6 indicate the presence of a strong 
cooling flow in the central r < 100 kpc region of the cluster. 
We have therefore also examined the spectral data for this 
region using a series of more sophisticated, multiphase mod- 
els in which the emission properties of the cooling flow were 
explicitly accounted for. The first such model (model CI) in- 
troduced an extra emission component into model A, with 
a spectrum appropriate for gas cooling at constant pressure 
from the ambient cluster temperature (following the pre- 
scription of Johnstone et al. 1992). The normalization of 
this component was parameterized in terms of a mass depo- 
sition rate, M, which was a free parameter in the fits. In the 
second case, model C2, the cooling gas was modelled as an 
isothermal cooling flow, following Nulsen (1998; we assume 
a value for rj — 1, where the integrated mass deposition rate 
within radius r, M oc r v ). The mean gas temperature and 
metallicity in both the cooling flow and isothermal emis- 
sion component (which accounts for the emission from gas 
at larger radii viewed in projection) were assumed to be 
equal. Finally, we also examined a more general emission 
model, model D, in which the cooling gas was modelled by 
a second, cooler isothermal emission component, with the 
temperature and normalization of this component included 



as free fit parameters. Model D provides a more flexible pa- 
rameterization, with an additional degree of freedom over 
models CI and C2, and invariably provides a good match to 
the more specific cooling-flow models at the spectral reso- 
lution and signal-to-noise ratios typical of ACIS cluster ob- 
servations. However, the parameter values determined with 
model D were not well constrained for Abell 2390, and thus 
we do not quote explicit results for this model here. 

With each of the cooling-flow emission models, we have 
also examined the effects of including extra absorption, us- 
ing a variety of different absorption models. In the first case 
(absorption model i) , the only absorption included was that 
due to cold gas in our Galaxy, with the equivalent column 
density fixed to the nominal Galactic value (Dickey & Lock- 
man 1990; For a single temperature emission model, this is 
identical to spectral model A). In the second model (model 
ii), the absorption was again assumed to be due to Galactic 
(zero redshift) cold gas, but with the column density, A^h, 
included as a free parameter in the fits. (For a single temper- 
ature emission model, this would be equivalent to spectral 
model B.) In the third case (absorption model iii), an intrin- 
sic absorption component with column density, AAh, due 
to cold gas at the redshift of the cluster was introduced. The 
absorber was assumed to lie in a uniform screen in front of 
the cooling flow, with the column density included a free 
fit parameter. In the fourth case (model iv), the intrinsic 
absorption was assumed to cover only a fraction, /, of the 
emission from the cooling flow. The fifth and final absorp- 
tion model (model v) was similar to model (iii) but with the 
gaseous absorber replaced by an intrinsic absorption edge, 
with the edge depth, r, and energy, E e d se , free parameters in 
the fits. This more general absorption model may be used to 
approximate the effects of a dusty and/or ionized absorber. 

In those cases where the absorption has been quantified 
in terms of an equivalent hydrogen column density, solar 
metallicity in the absorbing gas is assumed. We note that 
in absorption models (iii-v), the absorption acting on the 
ambient cluster emission was fixed at the nominal Galactic 
value. However, allowing the Galactic absorption to vary 
from this value did not significantly improve the fits. 

4.3 Results from the single-phase analysis 

The best-fit parameter values and la (A\ 2 = 1.0) and 90 
per cent (A^ 2 = 2.71) confidence limits determined from 
the fits in the 0.5 — 7.0 keV band with the single temper- 
ature models are summarized in Tables 1 and 2. The tem- 
perature profile determined with spectral model A is shown 
in Fig. ^|(a). The measured temperature is approximately 
isothermal beyond a radius of 200 kpc, out to the limits of 
the data at r ~ 1 Mpc. A combined fit to the data in the 
0.2 — 1.0 Mpc range with model A gives a mean tempera- 
ture of 11.5+^3 keV. Fitting the 0.2 - 1.0 Mpc results with 
a simple power-law model of the form kT oc r a we measure 
a = 0.0±0.2 (lcr bootstrap errors obtained using the Akritas 
& Bershady 1996 modification of the ordinary least squares 
statistic.) We observe a clear drop in the emission weighted 
temperature within the central 100 kpc, with a value for the 
central r = 50 kpc of 5.58t° 5 4 l keV. 

We detect marginal evidence for a metallicity gradient 
in the cluster, with a mean, best-fit value for the central 100 
kpc of Z = O.48+°;ioZ , which compares to a mean value 
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of 0.23±£;i3% for the outer 0.2 - 1.0 Mpc region (Fig. |a; 
la errors). The results determined with spectral model B in 
the 0.5 — 7.0 keV band also provide marginal evidence for 
increased absorption towards the cluster core (Fig. ^b). 

In all cases, the x 2 values obtained from the fits to the 
0.5 — 7.0 keV data with the single-temperature models are 
acceptable, the only marginal case being the data for the 
central 50 kpc region, where the spectrum is complicated by 
the effects of the cooling flow. 

4.4 Multiphase analysis of the cooling core 

The results determined from the more detailed, multiphase 
analysis of the central 100 kpc radius in the extended 
0.3—7.0 keV energy band are summarized in Table 3. The re- 
sults demonstrate a clear requirement for excess absorption 
in this region, over and above the nominal Galactic value 
(Dickey & Lockman 1990), using each of the different emis- 
sion models. The systematic variations between the column 
density measurements obtained with the different emission 
and absorption models are similar to those determined from 
previous ASCA studies {e.g. Allen 2000). Unfortunately, the 
present Chandra data for Abell 2390 cannot statistically dis- 
criminate between the single-phase and multiphase emission 
models for the central 100 kpc region, which provide com- 
parable % 2 values with each absorption model. (We note, 
however, that at some level the spectrum for the central 100 
kpc must be multiphase simply due projection effects, given 
the results on the temperature profile shown in Fig. |H|(a).) 

For our preferred cooling-flow emission models includ- 
ing intrinsic absorption, which provide the most consistent 
physical description of the spectral and imaging data for 
the central 100 kpc region (Section 6), we measure an in- 
tegrated mass deposition rate M = 200 — 300 M yr _1 . 
Fig. |(b) shows the projected temperature profile for Abell 
2390, corrected for the effects of the cooling flow using spec- 
tral model C2(iii). (For r > 100 kpc the results determined 
with spectral model A in the 0.5 — 7.0 keV band have been 
used). We see that correcting for the effects of the cooling 
flow does not have a major impact on the temperature pro- 
file measured in the central regions of the cluster. Assuming 
that the isothermal cooling flow model provides a reasonable 
description of the data, the drop in the central temperature 
shown in Figs. |^(a,b) should then reflect the mass distribu- 
tion in the cluster core. 

The data for Abell 2390 do not provide firm constraints 
on the nature of the intrinsic absorption in the cluster. Using 
the partial-covering absorption model (iv) , we find that high 
covering fractions are preferred. For the edge- like absorption 
model (v), the lower limit on the edge energy is essentially 
unconstrained: for emission models Cl(v) and C2(v), how- 
ever, the upper limit on the edge energy is inconsistent with 
the OIK edge of oxygen (£ c d gc ~ 0.54 keV), suggesting 
that the absorption is unlikely to be due to oxygen rich dust 
grains (e.g. Arnaud & Mushotzky 1998; Allen et al. 2000b) 
or a warm, ionized absorber (Buote 2000) . 

Assuming that the intrinsic absorption is due to cold gas 
lying in a uniform screen in front of the cooling flow, using 
spectral model C2(iii), we measure A7V H = 2.2+J;g x 10 21 
atom cm -2 . The mass of absorbing gas implied by this model 
is then M a b s ~ 3 x 10 7 r 2 bs AA r H Mq, where r a b s is the radial 
extent of the absorber in kpc and ATVh is the equivalent 



hydrogen column density in units of 10 21 atom cm 2 . For 
r abs ~ 50 kpc (Section 6.2) we obtain M abs ~ 1.7±S;| X 10 11 
M Q ( although Allen & Fabian 1997 and Wise & Sarazin 
2000 argue that for a geometry in which the absorbing ma- 
terial is distributed throughout the X-ray emitting region, 
the true mass is likely to be a few times higher). This mass 
is in reasonable agreement with the mass expected to have 
been accumulated by the cooling flow within the same radius 
over its lifetime; ~ Mt/2 ~ 2.8 ± 1.4 x 10 11 M (the factor 
two in the denominator arises from the assumption that the 
integrated mass deposition rate, M, grows approximately 
linearly with time). 

4.5 Spectral deprojection analysis 

The results discussed in the previous subsections are based 
on the analysis of projected spectra. In order to determine 
the effects of projection on the spectral results, we have also 
carried out a simple deprojection analysis of the Chandra 
spectral data. 

For this analysis we have used the same annular regions 
and have assumed that the emission from each spherical 
shell (the shells are defined by the same inner and outer 
radii as the annular regions) is isothermal and absorbed by 
the Galactic column density (spectral model A). The fit to 
the outermost annulus is used to determine the temperature 
and emission measure in the outermost spherical shell. The 
contribution from that shell to each inner annulus is then 
determined by purely geometric factors (e.g. Kriss, Cioffi & 
Canizares 1983). The fit to the second annulus inward is used 
to determine the parameters for the second spherical shell, 
and so forth, working inwards. Thus, the spectral model for 
the ith annulus working inwards is the weighted sum of i 
absorbed, isothermal models. This parallels the usual image 
deprojection procedure (e.g. Fabian et al. 1981). 

We have used the XSPEC code and Chandra data in the 
0.5 — 7.0 keV band. The data for all eight annular spectra 
were fitted simultaneously in order to correctly determine 
the parameter values and confidence limits. The spectral 
model used therefore has 2n + 1 free parameters (where n 
is the total number of annuli), corresponding to the tem- 
perature and emission measure in each spherical shell and 
the overall emission-weighted metallicity (the metallicity is 
linked to the same value at all radii) . Note that we have have 
not attempted to correct for residual emission from gas at 
radii beyond the outermost annulus since the steeply rising 
surface brightness profile of the cluster causes this emission 
to have a negligible affect on the results. 

The temperature profile determined with the spectral 
deprojection method is shown in Fig. [?]. 

4.6 Comparison with previous work 

The mean ambient temperature for Abell 2390 of kT — 
11.5^ 3 keV, determined from the combined analysis of 
the data in the 0.2 — 1 Mpc range, is in good agreement 
with the previous results of kT = 14.5±5 5 2 5 keV fr om Allen 
(1998) based on ASCA observations, and kT = ll.llJJ 
keV from Bohringer et al. (1998), based on a joint anal- 
ysis of ASCA and ROSAT Position Sensitive Proportional 
Counter (PSPC) data. (Both Allen 1998 and Bohringer et al. 
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1998 accounted for the effects of the central cooling flow in 
their modelling of the integrated cluster spectra. Bohringer 
et al. (1998) also examined a simpler, single-temperature 
emission model with which they measured a lower mean 
emission-weighted temperature of kT ~ 9 keV.) The am- 
bient temperature for Abell 2390 measured with Chandra is 
also in good agreement with the predicted value of ~ 12.0 
keV using the cooling- flow corrected feTx/I/Boi relation of 
Allen & Fabian (1998; we assume a bolometric luminosity 
of ~ 1.0 x 10 46 ergs -1 as measured by ASCA since the 
Chandra observations do not cover the whole of the clus- 
ter). 

The best-fit mass deposition rate from the cooling flow 
of M ~ 200 - 300 M Q yr" 1 , determined from the Chan- 
dra spectrum for the central 100 kpc region, is lower than 
previous measurements based on the analysis of integrated 
spectra for the whole cluster from ASCA (M = 1500t_f? 00 
Moyr" 1 ; Allen 2000), joint ASCA/ROSAT (M = 700±^° 
Moyr -1 ; Bohringer et al. 1998) and Beppo-SAX (M = 
700± 4 oo M yr" 1 ; Ettori, Allen & Fabian 2000) observa- 
tions (although the results are marginally consistent at the 
~ 95 per cent confidence level). In part, this difference is 
likely to be due to the fact that the ambient gas temper- 
ature in the centre of the cluster, corrected for the effects 
of the cooling flow, is lower than the mean value measured 
at larger radii (Fig. ^b). If this drop in the central ambient 
temperature is not accounted for (as was the case in the pre- 
vious ASCA and BeppoSAX studies, which could not spa- 
tially resolve the cooling flow from the hotter, outer cluster 
gas) then the cooler, ambient gas in the cluster core will also 
tend to be modelled as part of the cooling flow, and the total 
mass deposition rate will be overestimated. This effect is il- 
lustrated by the fact that a fit with spectral model C2(iii) to 
a single Chandra spectrum covering the entire central 500 
kpc (radius) of the cluster, gives M = 621+ 12 :!, Mq yr -1 , 
A Ah = l-6±g;| x 10 21 atom cm" 2 and kT = U.9±f% keV, in 
good agreement with the previous ASCA, ROSAT and Bep- 
poSAX results. However, a fit to the central 200 kpc radius 
region (the maximum possible size of any cooling flow; Sec- 
tion 6) gives M = 280 ± 110 M yr" 1 (kT = 8.6± | keV), 
in good agreement with the value listed in Table 3. These 
results, and the consistent findings from the spectral, image 
deprojection and X-ray colour profile analyses of the inner- 
most lOOkpc presented in Section 6, highlight the need for 
detailed spatially-resolved spectroscopy when attempting to 
study the properties of cooling flows in distant clusters. 

Finally, we note that the excess column density act- 
ing on the cooling flow component measured with spectral 
model Cl(iii) of A Ah = 2.4+,!,;;; x 10 21 atom cm" 2 is consis- 
tent with the previous measurement of 2.9+j'g x 10 21 atom 
cm - from ASCA, using the same model (Allen 2000). 



* Note that although the fit with spectral model C2(iii) to the 
Chandra spectrum for the central 500 kpc (radius) region overes- 
timates the mass deposition rate from the cooling flow, the mea- 
sured temperature is in good agreement with the true value at 
large radii in the cluster (Fig. pp). Thus, previous ASCA and 
Beppo-SAX studies which attempted to account for the effects 
of cooling flows on integrated cluster spectra (e.g. Allen 1998; 
Bohringer et al. 1998) may have overestimated the mass deposi- 
tion rates from the cooling flows but are likely to have provided 
reasonable estimates of the mean cluster temperatures. 



5 MEASUREMENT OF THE CLUSTER MASS 
PROFILE 

5.1 The mass model 

The observed X-ray surface brightness profile (Fig. |^a) and 
deprojected temperature profile (Fig. may together be 
used to determine the X-ray gas mass and total mass profiles 
in the cluster. For this analysis, we have used an updated 
version of the image deprojection code developed in Cam- 
bridge (see e.g. White, Jones & Forman 1997 for details). 
A variety of simple parameterizations for the cluster mass 
distribution were examined, to establish which could pro- 
vide an adequate description of the Chandra data. For those 
mass models providing reasonable fits, the best-fit param- 
eter values were determined using a simple iterative tech- 
nique.^ Spherical symmetry and hydrostatic equilibrium are 
assumed throughout. 

We find than a good fit (Xmm = 6.1 for 6 degrees of free- 
dom; hereafter DOF) to the Chandra data can be obtained 
using an NFW mass model 



p(r) 



Pcrit<5 c 



(r/r s ) (1 + r/r s 



where p(r) is the mass density, p c 
ical density for closure and 



(1) 



3Hy8nG is the crit- 



ic = 



200 



3 [ln(l + c)-c/(l + c)]' 



(2) 



with a scale radius, r s = 0.88_ ' 4g Mpc and a concentration 
parameter, c = 4.1+ 2 'y (68 per cent confidence limits). The 
normalization of the mass profile may also be expressed in 
terms of an equivalent velocity dispersion, a — \/50Hor B c 
(with r s in units of Mpc). The equivalent velocity dis- 
persion associated with the best-fit X-ray mass model of 
a = 1275+340 krns" 1 (68 per cent confidence limit), is in 
good agreement with the robust, optically-determined value 
of 1262lgg kms -1 (Borgani et al. 1999; from a re-analysis 
of the data of Carlberg et al. 1996 data). 

We note that the best-fit values of c, a and r a are 
correlated. Fixing r B — 0.88 Mpc in the NFW model, we 
determine a concentration parameter c = 4.10 ± 0.18 and 
an equivalent velocity dispersion a = 1275 ± 55kms _1 (68 
per cent confidence limit keeping r s fixed.) For r s = 0.4 
Mpc, we obtain c = 6.61 ± 0.28 and a = 935 ± 40 krns" 1 . 
Fixing r s — 2.0 Mpc, we obtain c = 2.53 ±0.11 and 
a = 1790 ± 80kms _1 . Note, however, that the mass dis- 
tributions within the central IMpc radius are similar in all 
three cases. 



t Given the observed surface brightness profile and a particular 
parameterized mass model, the deprojection code is used to pre- 
dict the temperature profile of the X-ray gas. (We use the median 
model temperature profile determined from 100 Monte-Carlo sim- 
ulations.) The predicted temperature profile is rebinned, using an 
appropriate flux weighting, and compared with the results from 
the spectral deprojection analysis (Fig. |?j). The x 2 difference be- 
tween the observed and predicted temperature profiles is then cal- 
culated. The parameters for the mass model are stepped through 
a regular grid of values in r s and a (or alternatively r B and c) to 
determine the best-fit values (which give the minimum x 2 ) and 
confidence limits. For Abell 2390 the final grid used ranges from 
0.1 - 2.5 Mpc in r s and 500 - 2500 krns -1 in a. 
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The best-fit NFW mass model, with r s = 0.88±J;|| Mpc 
and c = has a virial radius r2oo = cr s = 3.6 jli'o 

Mpc and an integrated mass within this radius, M200 = 
2.7lf;g x 10 15 M Q . The relationship between M 20 o, r s and 
c for Abell 2390 is consistent with that expected for such a 
massive cluster formed at a redshift z ~ 0.5 — 1.0 {e.g. Eke, 
Navarro & Frenk 1998) 

The deprojected X-ray gas temperature profile implied 
by the best-fitting NFW mass model (given the observed 
surface brightness profile) is shown overlaid on the depro- 
jected spectral results in Fig. ^ The model results have been 
rebinned to the same resolution as the spectral data. 

Finally, we note that the NFW mass model provides a 
significantly better-fit to the Chandra data for Abell 2390 
than a singular isothermal sphere (p(r) oc r~ 2 : Xmin = 33.1 
for 7 DOF). However, a variety of other, two-parameter 
models including a softened isothermal sphere (p(r) oc 
(1 + (r/r c ) 2 ): Xmin = 5.8 for 6 DOF with r c ~ 75 kpc), 
a King approximation to an isothermal sphere (p(r) oc 
(1 + (r/r c ) 2 )- 3/2 : Xmin = 6.2 for 6 DOF with r c ~ 150 
kpc), and a full isothermal sphere (Equation 4-125 of Bin- 
ney & Tremaine 1987: x min = 5.9 for 6 DOF with r c ~ 150 
kpc ) also provide good descriptions of the Abell 2390 mass 
profile. 

5.2 Comparison of X-ray and lensing mass 
measurements 

Abell 2390 is one of the best studied lensing clusters (e.g. 
Pello et al. 1991; Kassiola, Kovner & Blandford 1992; 
Narasimha & Chitre 1993; Pierre et al. 1996; Squires et al. 
1996; Bezecourt & Soucail 1997; Frye & Broadhurst 1998; 
Pello et al. 1999). The cluster exhibits an unusual, strongly 
lensed 'straight arc' approximately 38 arcsec (174 kpc) away 
from the nucleus of the central galaxy (Pello et al. 1991) in 
addition to many other arcs and arclets (e.g. Bezecourt & 
Soucail 1997; Pello et al. 1999). Pierre et al. (1996) present 
a two-component mass model for the central regions of the 
cluster and measure a projected mass within the radius de- 
fined by the brightest arc of 1.6 ±0.2 x 10 14 M©. Squires et 
al. (1996) present a weak lensing analysis of the cluster and 
determine an azimuthally averaged mass profile covering the 
central ~ 4.3 arcmin (1.2 Mpc). 

Fig. ^| shows the projected mass profile determined from 
the Chandra X-ray data, with the Squires et al. (1996) weak 
lensing results and Piere et al. (1996) strong lensing results 
overlaid. For the X-ray analysis, we assume that the NFW 
mass models extend out to V2oo in each case. (The lim- 
its on the X-ray results are the maximum and minimum 
masses at each radius for the range of NFW models with 
X 2 < Xmin + 2.30; the 68 per cent confidence contour in the 
r s — a plane). The agreement between the X-ray and lens- 
ing mass results in Fig. |^ is reasonable at all radii studied. 
The mean scatter of the lensing results about the best-fit 
X-ray mass profile within the central 1 Mpc region is < 20 
per cent. The agreement between the independent lensing 
and X-ray mass measurements, together with the consistent 
results on the equivalent X-ray and observed optical galaxy 
velocity dispersions, confirms the validity of the hydrostatic 
assumption used in the X-ray analysis and suggests that the 
mass profile in Abell 2390 has been robustly determined. 

We note that at small radii, the strong lensing mass 



of Pierre et al. (1996) slightly exceeds (at ~ la signifi- 
cance) the best-fit value determined from the Chandra data 
(1.19±[J-jf x 10 14 M© within r = 174 kpc). This difference, 
albeit marginal, may be related to the residual substructure 
on these scales seen in the X-ray (Section 3.1) and optical 
images and strong lensing mass map (Pierre et al. 1996), 
and could indicate a slight enhancement of the central lens- 
ing mass due to an alignment of the dominant mass clumps 
in the cluster and/or the presence of additional, non-thermal 
pressure support of the X-ray gas in the cluster core. 

5.3 The X-ray gas mass fraction 

The X-ray gas-to-total-mass ratio as a function of radius, 

tas(f), determined from the Chandra data is shown in Fig. 
We find that the best-fit / gas value rises with increasing 
radius with / gas = 0.21±%™hm' S at r = 0.9 Mpc (conser- 
vative limits determined by combining the la errors on the 
integrated gas mass at each radius with the uncertainties in 
the total mass distribution shown in Fig. ^ in quadrature). 
This value is consistent with the previous measurement of 
17.5 ± 2.4 per cent (la limits) at r = IMpc from Ettori 
& Fabian (1999) using ROSAT PSPC and ASCA data and 
assuming strict isothermality in the X-ray gas. 

Following the usual arguments, which assume that the 
properties of clusters provide a fair sample of those of 
the Universe as a whole (e.g. White et al. 1993; White 
& Fabian 1995; Evrard 1997; Ettori & Fabian 1999; Bah- 
call et al. 1999), we may use our result on the X-ray gas 
mass fraction in Abell 2390 to estimate the total matter 
density in the Universe, f2 m . Assuming that the luminous 
baryonic mass in galaxies in Abell 2390 is approximately 
one fifth of the X-ray gas mass (e.g. White et al. 1993; 
Fukugita, Hogan & Peebles 1998) and neglecting other pos- 
sible sources of baryonic dark matter in the cluster, we ob- 
tain il m = (fib/1. 2/g as ), where Qb is mean baryon density 
in the Universe and hso is the Hubble constant in units of 
SOkms -1 Mpc" 1 . For Qb/ifo = 0.0820±0.0072 (O'Meara et 
al. 2000), fi m = 0.33 ± 0.16^ Q 5 |] 

6 THE PROPERTIES OF THE COOLING 
FLOW 

6.1 Radial properties of the cluster gas 

The results on the electron density and cooling time as a 
function of radius, determined from the image deprojection 
analysis using the best-fit NFW mass model, are shown 
in Fig. |Io| Within the central 500 kpc radius, the elec- 
tron density profile can be parameterized (x 2 = 29.1 for 
52 degrees of freedom) by a /3-profile with a core radius, 
r c = 21 ± 3 kpc, /3 = 0.378 ± 0.014 and a central density, 
n e (0) = 8.5±0.8x 10" 2 cm" 3 (la errors). The core radius for 
the electron density distribution is slightly smaller than the 
value measured directly from the projected surface bright- 
ness profile under the assumption of strict isothermality in 

t Accounting for any additional, dark baryonic matter in the clus- 
ter would lower the measured value of f2 m . Likewise, if the value 
of /gas increases towards larger radii in Abell 2390, the true value 
for £7 m will be lower than our quoted result. 
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the X-ray gas (Section 3.2). Indeed, the evidence for any flat 
central core in the electron density distribution is marginal 
and a simple broken power-law model, with rbroak = 46_ 11 
kpc and slopes interior and exterior to the break radius of 
-0.48 ± 0.05 and -1.13 ± 0.04, provides as good a fit to 
the electron density profile in the central 500 kpc region 
(X 2 = 30.4 for 51 de grees of freedom). 

For an assumed Galactic column density of 6.8 x 
10 20 atom cm 2 , we measure a central cooling time (i.e. the 
mean cooling time within the central 1.97 arcsec or (~ 9 
kpc) bin of t CO oi = 3.9lJ' g x 10 8 yr, and a cooling radius, 
at which the cooling time first exceeds a Hubble time, of 
'"cool = 175lg° kpc. (Errors on the central cooling time are 
the 10 and 90 percentile values from 100 Monte Carlo simu- 
lations. The upper and lower confidence limits on the cooling 
radii are the points where the 10 and 90 percentile values 
exceed and become less than the Hubble time, respectively.) 

6.2 X-ray colour profile analysis 

We have constructed an X-ray 'colour' profile for the clus- 
ter in order to determine the size of the central region in 
which significant cooling occurs. Two separate images were 
created in the energy bands 0.5 — 1.3 and 1.3 — 7.0 keV 
(0.62 — 1.60 and 1.6 — 8.6 keV in the rest frame of the 
source) with a 1.97 arcsec (4 raw detector pixels) pixel scale. 
These soft and hard X-ray images were background sub- 
tracted and flat fielded (taking full account of the spectral 
energy distributions of the detected photons). All significant 
point sources were masked out and excluded from the analy- 
sis. Azimuthally- averaged surface-brightness profiles for the 
cluster were then constructed in each energy band, centered 
on the overall peak of the X-ray emission (Section 3). The 
X-ray 'colour' profile formed from the ratio of the surface 
brightness profiles in the soft and hard bands is shown in 

From examination of Fig. hll we see that at large radii 
the observed X-ray colour ratio is approximately constant, 
with a mean value of 1.01 ±0.03 (la error determined from 
a fit to the data between radii of 150 — 300 kpc). By com- 
parison with simulated spectra we find that this result is 
consistent with an isothermal plasma with a temperature 
kT = 11.0 ± 1.5 keV. (We assume a metallicity Z ~ OAZq 
and Galactic absorption.) Within a 'break' radius of 70t 2 2 
kpc (la errors determined from a % 2 fit with a broken power- 
law model), however, the colour ratio rises sharply, indicat- 
ing the presence of significantly cooler gas. 

6.3 Analysis of the mass deposition profile 

The outermost radius at which cooling occurs may also be 
expected to be associated with a 'break' in the X-ray surface 
brightness profile and, more evidently, the mass deposition 
profile determined from the deprojection code. The mass 
deposition profile from the cooling flow, which is a parame- 
terization of the X-ray luminosity distribution in the cluster 
core (see e.g. White et al. 1997), is shown in Fig. ^| Fit- 
ting this profile with simple a broken power-law model, we 
determine a break radius of 51 ± 12 kpc (la errors), in rea- 
sonable agreement with the break radius determined from 
the fit to the X-ray colour profile (Section 6.2) and the elec- 
tron density distribution (Section 6.1) . The best-fit broken 



power-law model is shown overlaid on the mass deposition 
profile in Fig. |l2|. 

The slopes of the mass deposition profile, internal and 
external to the break radius are 1.72 ± 0.16 and 0.76 ± 0.16, 
respectively. Accounting only for absorption due to cold 
gas with the nominal Galactic column density, we deter- 
mine an integrated mass deposition rate within the break 
radius of 212+^ M Q yr -1 . If we also account for the pres- 
ence of intrinsic absorption, with the properties determined 
using spectral model C2(iii), the mass deposition rate within 
the break radius rises to 247+gi Moyr" 1 , in good agree- 
ment with the spectral result for the central 100 kpc of 
2421" M Q yr" 1 . 



6.4 The age of the cooling flow 

Allen et al. (2000a) discuss a number of methods which may 
be used to estimate the ages of cooling flows from X-ray 
data. (Such ages are likely to relate to the time intervals 
since the central regions of clusters were last disrupted by 
major subcluster merger events.) Essentially, these methods 
identify the age of a cooling flow with the cooling time of 
the X-ray gas at the break radius in either the X-ray colour 
or deprojected mass deposition profile. 

Using Fig. [lo](b) , we see that cooling time of the X-ray 
gas at the break radius in the X-ray colour profile (70t 2 2 
kpc) lies in the range 3.21q j? Gyr. (The cooling time at the 
break radius is measured from a least-squares fit to the data 
in Fig. [lo](b) over the 10 — 100 kpc range using a power- 
law model.) If we instead identify the age of the cooling 
flow with the cooling time at the break radius in the mass 
deposition profile in Fig. [12l (51 ± 12 kpc), we infer an age 
for the cooling flow of 2.3Iq^ Gyr. (In both cases we assume 
that no intrinsic absorption acts beyond the outer edge of 
the cooling flow, which is reasonable if the absorbing matter 
is accumulated by the flow.) 

In summary, we see the X-ray colour profile, image de- 
projection analysis and spectral data provide consistent re- 
sults on the properties of the cooling flow in Abell 2390, 
indicating a mass deposition rate in the range 200 — 300 
Mq yr -1 and an age of 2 — 3 Gyr. 

6.5 On the rising ambient temperature profile 
within the cluster core 

In principle, the results on the X-ray gas temperature and 
surface brightness within the cluster core may be used to 
distinguish between an NFW model for the dark matter 
distribution and alternative models with a steeper central 
cusp (e.g. Moore et al. 1999). The results for Abell 2390 
presented here are certainly consistent with an NFW profile 
(see Fig. Q). However, given the complexity of the gas within 
the central 100 kpc and the relatively short exposure time 
of the present observation (which limits the number of inde- 
pendent spectra which can be extracted from the region of 
the cluster core), we do not attempt to explore such models 
here. This issue will be better addressed with Chandra data 
for nearer, brighter systems. 

The fact that the ambient temperature profile, cor- 
rected for the effects of the cooling flow, rises with increasing 
radius throughout the central ~ 200 kpc, is interesting. In 
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detail, the density and temperature structure in a cluster 
core will depend on the thermal history of the gas, as well 
as the underlying dark matter distribution and, subject to 
pressure equilibrium and convective stability, may be flat, 
increase or decrease with radius. As discussed in Section 
4.6, the presence of relatively cool, ambient gas in the cen- 
tral regions of Abell 2390 may have lead to overestimates of 
the mass deposition rate in previous studies with ASCA and 
Beppo-SAX in which this region was not spatially resolved. 

There are several reasons why relatively cool, dense gas 
might exist beyond the outer edge of the present-day cool- 
ing flow in Abell 2390. The first is that a pre-existing cool- 
ing flow may have been disturbed several Gyr ago. Cooler, 
denser gas could then have been spread out over several 100 
kpc, with the cooling flow having so far only re-established 
itself within the inner ~ 50 — 70 kpc. A second possibility 
is that the densest gas from the cores of infalling subclus- 
ters may have been stripped and deposited over time in the 
core of the main cluster without strong shocking (Fabian & 
Daines 1991). The situation seen in the Chandra data for 
Abell 2142 (Markevitch et al. 2000) is reminiscent of this, 
with a sharp drop in temperature from kT ~ 14 keV to 
kT ~ 7 keV (and a corresponding rise in density) observed 
at r ~ 370 kpc, moving inwards from northwest of the clus- 
ter centre. Such a sharp drop in temperature at these radii 
cannot be due to any current cooling flow in the cluster. 



7 CONCLUSIONS 

The main conclusions from this work may be summarized 
as follows: 

(i) We have measured the distribution of mass in Abell 
2390 using Chandra X-ray observations. The mass profile 
can be well-modelled by an NFW profile with a scale ra- 



dius r B 



0.88^ 



Mpc and a concentration parameter, c = 



4.lli 7 (^S P er cen t confidence limits). The normalization of 
the mass profile may also be expressed in terms of an equiva- 
lent velocity dispersion, a = ^/50Hor B c = 1275^3^0 kms -1 , 
in good agreement with the optically-determined value of 
1262+^ kms- 1 (Borgani et al. 1999). 

(ii) The best-fit Chandra mass model is in good agree- 
ment with independent measurements from strong and weak 
lensing studies. The mean scatter between the X-ray and 
lensing values within the central r = 1 Mpc radius is < 20 
per cent. 

(iii) The X-ray gas to total mass ratio rises with increas- 
ing radius within the central 0.9 Mpc, with / gas = O^lj^g'^Q 
at r — 0.9 Mpc. Following the usual arguments, this re- 
sult on the X-ray gas mass fraction may be converted into 
a constraint on the mean mass density in the Universe, 
On = 0.33±0.16/i^ ' 5 . 

(iv) The X-ray gas temperature rises with increasing 
radius within the central r ~ 200 kpc, and then remains 
approximately isothermal with kT — 11.5ii'f keV out to 
r ~ 1.0 Mpc. 

(v) The azimuthally-averaged 0.3 — 7.0 keV surface 
brightness profile exhibits a small core radius and a clear 
break at r ~ 0.5 Mpc, where the slope changes abruptly 
from Sx °S r" 1 ' 5 to Sx °S r~ 3 6 . 

(vi) The best-fit mass deposition rate from the cool- 
ing flow, determined in a consistent manner from the spec- 



tral and imaging Chandra data, lies in the range 200 — 300 
Mq yr _1 . This value is lower than previous estimates based 
on integrated ASCA and Beppo-SAX spectra for the entire 
cluster, which could not resolve the drop in the central, am- 
bient gas temperature. We estimate an age for the cooling 
flow of 2 — 3 Gyr. 
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Table 1. The results from the analysis of the annular spectra using spectral model A in the 0.5 — 7.0 keV band. Temperatures (kT) are in 
keV, metallicities (Z) in solar units and normalizations (K) in units of (10 -17 /(47t(Da(1 + z)) 2 )) J n c rijidV, where D\ is the angular size 
distance to the source in cm, n c is the electron density (in cm" 3 ), and rijj is the hydrogen ion density (in cm -3 ). The absorbing column 
density has been fixed at the nominal Galactic value of 6.8 X 10 20 atom cm -2 . The total \ 2 values and number of degrees of freedom 
(DOF) in the fits are listed in column 5. Error bars at both the la (A\ 2 = 1.0) and 90 per cent (A% 2 = 2.71; in parentheses) confidence 
limits on a single interesting parameter are listed. For the outermost 0.7— 1.0 Mpc annulus, where the instrumental background is most 
significant, only data in the energy range 0.5 — 6.0 keV range were used. 
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Table 2. The results from the analysis of the annular spectra in the 0.5 — 7.0 keV band using spectral model B. The absorbing column 
density (Nr) is in units of 10 20 atom cm -2 . Other details as in Table 1. 



Model B kT Z N H X /DOF 

n 50 4 --+0.47(4-0.82) „ .„+0. 16(4-0.26) q „ fi +1.25(+2.08) 010/70 

- 50 4 - 77 -0.40(-0.63) °- 48 -0.14(-0.23) 9 - 86 -1.24(-2.00) 91.8/78 

n mn fi 1s +0.49(+0.83 „ „.+0.09(+0.16) o 71 +0.95(4-1.58) 1900/199 
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Table 3. The best-fit parameter values and l<x (and 90 per cent; in parentheses) confidence limits determined from the multiphase 
analysis of the central 100 kpc region in the extended 0.3 — 7.0 kcV band. The mass deposition rates from the cooling flow (M) are in 
units of Mq yr~ x . Column densities (JVh) and intrinsic column densities (AJVh) are in units of 10 20 atom cm -2 . The limits on the edge 
energy (E e dge) determined with absorption model (v) are in keV. Other details as in Table 1. 
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Figure 1. (Left panel) The raw 0.3 — 7.0 keV Chandra image of Abell 2390. The pixel size is 4 detector pixels (1.97 arcsec). (Right 
panel) Contour plot of the same region, adaptivcly smoothed using the code of Ebeling, White & Rangarajan (2000), with a threshold 
value of 3.5o\ The contours have equal logarithmic spacing. 
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Figure 2. An adaptivcly smoothed 0.3 — 7.0 keV image of the 
central regions of Abell 2390 on a finer spatial scale (pixel size 
0.492 X 0.492 arcsec 2 , which is equivalent to 1 X 1 raw detector 
pixels). Note the ridge of enhanced emission extending ~ 3 — 4 
arcsec to the northwest of the X-ray peak which is coincident with 
the excess blue continuum and optical/UV line emission reported 
by Lemonon et al. (1998) and Hutchings & Balogh (2000). 
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Figure 3. (a) The background-subtracted, fiat-fielded, azimuthally-averaged radial surface brightness profile for Abell 2390 in the 0.3—7.0 
keV band. The binsize is 0.984 arcsec (4.65 kpc). (b) The same profile, rebinned by a factor of 4, in the range 80 — 900 kpc, with the 
best-fitting broken power-law model overlaid. At the break radius of 491^^. 6 kpc, the slope of the surface brightness profile (£?x (x r a ) 
changes from a = -1.46 ± 0.02 to a = -3.60 ± 0.17. 
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Figure 4. The 0.3 — 7.0 keV surface brightness profile in quad- 
rants covering position angles 315 — 45 degrees (filled circles), 
45 — 135 degrees (open circles), 135 — 225 degrees (triangles) and 
225 — 315 degrees (squares). The extension towards the East is 
also evident in Fig. []]. The results for the other quadrants are in 
good agreement with each other. 
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Figure 5. (a) The projected X-ray gas temperature profile (and la errors) measured with spectral model A in the 0.5 — 7.0 keV energy 
range. The dotted lines mark the 90 per cent confidence limits on the mean ambient cluster temperature determined from a joint analysis 
of the data in the 0.2 — 1.0 Mpc region, (b) The ambient X-ray gas temperature as a function of radius, corrected for the effects of the 
cooling flow in the central lOOkpc using spectral model C2(iii; see Section 4.4). For r > 100 kpc, the values are determined in an identical 
manner to those in (a). 



© 0000 RAS, MNRAS 000, 000-000 



18 S. W. Allen et al. 




200 400 600 800 1000 
Radius (kpc) 



cv? 

I 



B 

V 



O 
■+-> 

(C 

o 

o 




200 400 600 800 1000 
Radius (kpc) 



Figure 6. (a) The (projected) variation of mctallicity in the cluster measured with spectral model A in the 0.5 — 7.0 keV energy range, 
(b) The absorbing column density measured with spectral model B. The dashed line shows the Galactic column density determined from 
HI studies (Dickey & Lockman 1990). 
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Figure 7. The predicted deprojected temperature profile (grey 
curve) determined from 100 Monte-Carlo simulations using the 
best-fitting NFW mass model (with r s = 0.88 Mpc, c = 4.1 and 
a = 1275 km s -1 ; see Section 5). The predicted profile has been 
binned to the same spatial resolution as the spectral deprojection 
results (solid points; Section 4.3) which are also shown overlaid. 
The agreement between the deprojected spectral results and best- 
fit NFW mass model predictions (reduced xi = 1-0 f° r v — 6 
degrees of freedom) indicates that the NFW mass model provides 
a good description of the spatially-resolved Chandra spectra. 
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Radius (kpc) 

Figure 8. A comparison of the projected total mass determined from the Chandra X-ray data (Section 5) with the strong lensing result 
of Pierre et al. (1996; filled circle) and the weak lensing results of Squires et al. (1996; open triangles). The best-fit NFW X-ray mass 
model is shown as the dotted curve (parameters r s = 880 kpc, a = 1275 km s -1 ). The limits on the X-ray results (the maximum and 
minimum masses at each radius for the range of NFW models with x 2 < xL in + 2.30; the 68 per cent confidence contour in the r s — a 
plane) are shown as solid curves. 
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Figure 9. The ratio of the X-ray gas mass to total gravitating 
mass as a function of radius. The three curves show the best- 
fit value (dotted curve) and conservative confidence limits (solid 
curves; determined by combining the la errors on the integrated 
gas mass at each radius with the uncertainties in the total mass 
distribution shown in Fig. ph. At r = 0.9 Mpc we measure / gas = 
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Figure 10. The results on (a) the electron density and (b) the cooling time, determined from the X-ray image dcprojcction analysis 
using the best-fit NFW mass model. Error bars are the la errors determined from 1000 Monte Carlo simulations. A Galactic column 
density of 6.8 X 10 20 atom cm -2 and a metallicity of 0.4 solar are assumed. 



© 0000 RAS, MNRAS 000, 000-000 



Chandra measurements of the distribution of mass in Abell 2390 23 



> 

CO 

M 
o 

CO 



> 

CO 

M 

CO 



I 

LO 

d 



LO 



cvj 



LO 



1 1 1 1 1 

= <- 

I 

-f H 


- 








100 200 
Radius (kpc) 



300 



Figure 11. The X-ray colour profile, formed from the ratio of the 
counts in the 0.5 — 1.3 and 1.3 — 7.0 keV bands, as a function of ra- 
dius. At large radii the ratio is approximately constant and consis- 
tent with a cluster temperature of 11. Oil. 5 keV. Within a 'break' 
radius of 70^ 2 j kpc, however, the colour ratio rises sharply. The 
grey curve shows the predicted colour profile for a cooling flow of 
age ~ 3 Gyr from the image deprojection analysis (Section 6). We 
use the best-fit NFW mass model and assume that the cooling gas 
is intrinsically absorbed by a column density of 2.2 X 10 21 atom 
cm -2 , as determined with spectral model Cl(iii). 
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Figure 12. The mass deposition rate (MDR) in units of 
Mq yr _1 , determined from the image deprojection analysis using 
the best-fit NFW mass model. The dotted line is the best-fitting 
broken power-law model (Section 6.3). Error bars are the 10 and 
90 percentile values from 1000 Monte Carlo simulations. A Galac- 



tic column density of 6.8 X 10 20 
solar arc assumed. 



atom cm , a metallicity of 0.4 
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